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SUMMARY 


MISS' 


The spectrum of a phase modulated carrier is analyzed. The 
analysis is based on the RMS phase deviation only, so that knowl- 
edge of the frequency spectrum of the modulating time function is 
not required. Upper and lower limits of the amplitudes of the spec- 
tral components of the phase modulated carrier are derived at, the 
main emphasis being given to the component of carrier frequency. 

The usefulness and simplicity of the equations derived at is demon- 
strated in the analysis of a carrier phase modulated by a pseudo 
random code and two phase modulated subcarriers. The calcula- 
tions are carried out without introducing approximations such as . ff 
substituting voice or other information by a single frequency. accaub 
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APOLLO USB-SYSTEM 

Multisignal Phase Modulation Problems and Solutions 


INTRODUCTION 

For the Apollo missions, both tracking and communication will be handled 
by the Unified S-band System (USBS). In this system, one main carrier is phase 
modulated by a base band consisting of a pseudo random range code and a number 
of phase modulated subcarriers [l]. In this paper the frequency spectrum of the 
modulated carriers is analyzed. Special attention is given to the calculation of 
the main carrier component because the success of the USBS hinges on the capa- 
bility of the receivers to lock on to and track the main carrier component. To 
analyze the power in the main carrier component under worst case conditions is 
therefore of vital importance. Although the analysis presented here was made 
with the USBS in mind, the results are, of course, applicable to any system utiliz- 
ing multisignal phase modulation. 

The conventional way of analysis of multisignal phase modulation is to expand 
the base band into a Fourier series and then to treat the problem as multi-tone 
phase modulation. A number of papers (see Reference [3], [4], [5], and [6]) have 
been written about multi-tone phase modulation. Unfortunately, the equations 
derived at are long and hard to evaluate as soon as the modulating time function 
consists of more than two or three frequency components. Another drawback of 
this method is that it is often difficult or impossible to expand the base band into a 
Fourier series. This is the case for a pseudo random code or a phase modulated 
subcarrier, which is modulated with an unknown periodic time function. However, 
the power of a pseudo random code is constant and the power of a phase modulated 
subcarrier is constant and independent of the modulating time function.* The 
power in the base band is thus easily calculated and is independent of the infor- 
mation modulated onto the subcarriers, A method of analysis of phase modula- 
tion, based on the power of the modulating time function (= base band) rather on 
the Fourier spectrum, has therefore been developed and is presented in this 
paper. With this type of analysis the exact value of the main carrier component 
cannot be calculated, but upper and lower limits can be calculated. This is not a 
limitation, but an advantage, because the upper and lower limits, rather than any 
particular value in between, are essential for circuit design and system analysis. 


This is only true if the Fourier spectrum of the modulating time function does not contain a 
component of the carrier frequency. 
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The analysis has been based on complex Fourier analysis because of the 
simplicity in derivations and lucidity of results obtained with the complex method. 
A brief recapitulation of complex analysis of phase modulation* is given in the 
next chapter, which may be omitted by the reader completely familiar with the 
subject. 

Only phase modulation has been considered in this paper because frequency 
modulation with f(t) is identical to phase modulation with f(t) dt. The 
integral of a function is more well behaved than the function itself and some time 
functions can therefore only be treated with phase modulation and not with fre- 
quency modulation. An example is the frequently used phase modulation with a 
square wave, also known as phase shift modulation. The concept of phase modu- 
lation thus completely covers frequency modulation. 


*The complex analysis of phase modulation is based on lectures given by Professor E. Lofgren in 
1952 at the Royal Institute of Technology, Stockholm, Sweden. 
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2. General Analysis of Multi-Tone Phase Modulation 




In this chapter we will analyze the case where a carrier cos (a) 0 t + \p 0 ) is 
phase modulated by a time function f(t). We assume that f(t) can be expanded 
in a Fourier series 


f(t) = 


n=N 



q> a sin (<“ B t + 


(2-1) 


where N may go to infinity. The D-C or average term of f(t) is included in the 
phase angle 0 O of the carrier and f(t) may thus be considered to be an A-C 
signal. If we normalize the peak amplitude of the modulated carrier to 1, then 
the instantaneous amplitude of the modulated carrier is 


V = cos (co Q t + 4> 0 + f(t)) 


( 2 - 2 ) 


or 


V = cos 


w 0 t + ^0 + 


n*N 



<p„ sin +'/'„) 


Using complex notation, we obtain 


V = 


j j"° 


t+ '/'0 + ^ < P „ sin (“„ t + '/ / n > 


»v 

(2-3) 


where 


v = Re{V} 


(2-4) 


In appendix A it is shown that V can be expanded in an infinite series 
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j(<v + 0O> 


jm(6> n t +l// n ) 


(2-5) 


V = e 


TT 


m=co 


z: 


JM e 


where J m are the Bessel functions of the first kind and order m . This is a very 
general expression for the frequency spectrum resulting from a phase modulation. 
Unfortunately, it is very difficult to obtain any information about the amplitudes 
of the frequency spectrum from Equ. (2-5) if N is a large number. In the following 
paragraphs, Equ. (2-5) will be analyzed for some cases of practical importance 
and simplified and useful expressions will be derived at. Before going into this, 
let us apply Equ. (2-5) to two special cases: 


Case A Let f(t) be a sinusoidal function 


f(t)=0sin(pt + y ) 


(2-6) 


For this case Equ. (2-5) reduces to 


V = 


,*<"0* + V'o> 


L 


m=-co 


T (>?) e im ( pt 


Using the relation [7] 


(2-7) 


= <-»>■ J„<x) 


we can rewrite Equ. (2-7) 


V = e 


i(<V +l /V 


J 0 (& + 2 J m (^ cos m (P t + ^ 
I n=2 e 


(2-8) 


uu 

* z: 

m*l 


+ i 2 / J m (<^>) sin m(pt + <//) 


4 




where 2 denotes summarization over even m and 2 n denotes summarization 
over odd m. Equation (2-8) can be interpreted by a phasor diagram utilizing 
the complex plain as shown in Fig. 1. The parenthesis in Equ. (2-8) is represented 
by the complex quanity Z , which is composed of infinitely many real and imagi- 
nary components. Z is of magnitude 1 due to the normalization introduced in the 
beginning of this chapter and performs an angular oscillation within the limits 
±$>. The factor e + '^o ) j s represented by the rotating time line t. In Equ. 
(2-2), we assumed a cosine function for v and v is therefore the projection of 
Z on the time line t . If we instead had assumed a sine function in Equ. (2-2) 
then v would be the projection of Z on a line perpendicular to the time line t . 


/ 



Figure 1-Phosor diagram for phase modulation of a carrier with one sinewave of 
angular frequency p and amplitude 4 > 
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Case B. Let f(t) be a subcarrier co x phase modulated with sin pt, thus 


A 


f(t) = <£ 0 sin (Wjt +\p t + ^ sin pt) 


(2-9) 


or 


f(t) = Im 


j(a > 1 +\p x + s in pt) 


> 


(2-10) 


The series'expansion of Equ. (2-10) is the same as Equ. (2-7) and we therefore 
have 


n*<P 

f(t) = $ 0 J n ) sin ((Wj + np) t + \p x ) (2-11) 

n—oo 


Equ. (2-11) is of the same form as Equ. (2-1) with 


=^0 J„(^) 


"n = + n P 


*n=*. 

Substitution into Equ. (2-5) yields 


V = 


H"o* +, /'o ) 


n 3 ® m «0D 


n*«H30 **-oo 


( 2 - 12 ) 


where 

4> 0 ~ peak phase swing of main carrier modulation 
4> x * peak phase swing of subcarrier modulation 
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coq « angular main carrier frequency 
^ = angular subcarrier frequency 
P = angular frequency of subcarrier modulation 
= constant phase angle 

This equation gives the complete frequency spectrum of a carrier, which is phase 
modulated with a phase modulated subcarrier. 

The complexity of Equ. (2-5) and (2-12) indicate that a better method of 
analysis is desirable, hi the next chapters the analysis will be based on the 
RMS value or average power 4 > rather than on the Fourier spectrum of the 
modulating time function. 


3. Phase Modulation with Harmonically Unrelated Frequencies 


If the multiplications of the sums in Equ. (2-5) are carried out, an expression 
of the form 


±A 0 

i=— oo 

is obtained, where 


< 3_2) 

n = l 


co. - ± M , co ± ± . . . ±M co ± . . . ±M n o; (3 3) 

i 1122 nn INN 

co = angular frequencies of the modulating function f ( t ) 
M n = an integer, including 0 

If w. is always different from zero if any one or several of the integers M n are 
different from zero, then the angular frequencies are harmonically unrelated. 

We see from Equ. (3-1) that in this case the amplitude of the component of car- 
rier frequency is A 0 . From Equ. (2-5) we see that co. is zero only for m ~ 0 
(provided, of course, that co + np ¥ 0, see footnote on page 1). The carrier 
component amplitude can thus be obtained from Equ. (2-5) by putting m = 0 , 
which leads to Equ. (3-2). An example of a time function with harmonically 
unrelated frequency spectrum is a subcarrier phase -modulated by 4> x sin pt . 

The frequency spectrum for this time function was given by Equ. (2-12). The 
carrier component is obtained by putting m = 0 in Equ. (2-12). The resulting 
expression can be brought in the same form as Equ. (3-2) after change of 
variables. 
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If, on the other hand, the frequencies are harmonically related, the amplitude 
of the carrier component is no longer given by Equ. (3-2). This case will be 
analyzed in the next chapter. 

If the power of the modulating time function (= base band) is $ , then 


N 



n = l 


(3-4) 


Under this condition it is shown in Appendix B, Equ. (B-12) and (B-15), that 


n«N 

J„ C/2 0) £ j (" J 0 (tj5„) — 

ns] 



(3-5) 


provided all cp n 1 5.1 for N > 2. If N = 1 then equ. (3-5) reduces to an identity, 
which is true for all values of <p n . It can further be shown, (see Appendix B) that 


lim 

N-oo 



= e -* 2 /2 


(3-6) 


J 0 (/2 <£)is always positive in the region 0 £ /2 — 2.4048, where 2.4048 is the 

first zero of J 0 (/2 <£)> and we can combine Equ. (3-2), (3-5), and (3-6) into 


J 0 (/2 4>)< A 0 i 


>- 4> 2 /2 


0±V2<P Z 2.4048 


(3-7) 


For /2 4> > 2.4048 we obtain instead 


0 1 A 


< 
o - 


riJo(^2 4 >)\ 

1 e -^ /2 




whichever is greater 




(3-8) 


2.4048 5 /2 0 i 5.1 /2 


The upper limit for <p is obtained from equ. (3-4) for N = 2 and tp n = 5.1. 

Equations (3-7) and (3-8) thus show that upper and lower limits of the carrier 
component can be calculated if only the power of the modulating time function is 
known. The lower limit is obtained if the modulating power is contained in one 
modulating frequency and the upper limit is obtained if the power is divided 
equally among infinitely many modulating frequencies. 

The Equ. (3-7) and (3-8) are shown graphically in Fig. 2. It is seen that for 
<p< 1 the upper and lower limits are very close together. The amplitude of the 
carrier component varies therefore only very little in this region if the frequency 
spectrum of the modulating time function is changed. 

The reduction of the main carrier component due to modulating is commonly 
referred to as modulation loss. The modulation loss L for the carrier can be 

m 

defined as the ratio of total output power and the power of the carrier component 
Expressing L m in db we obtain 


L„ = 10 log 


Total output power 
Carrier component power 


db 


Using the limits for the carrier component given in Equ. (3-7) and (3-8) we obtain 


10 log 1 > L > 10 log e^ (3-7 a) 

J 0 (^ 4>) 2 - m - 
0 <^2 4 >< 2.4048 
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RADIANS 


Figure 2-Upper limit Aq ma x an< J lower limit Aq m j n of the carrier component amplitude Aq. 
The carrier is phase modulated with RMS phase swing 0 by a time function f(t), whose 
frequency spectrum is not harmonically related. 


and 


L m > J 


10 log 


10 log 


J 0 (^2 0)2 


whichever is smaller 


(3-8a) 


for 


2.4048 < V2 4 > < 5.1 /2 

Equ. (3-7a) is shown graphically in Fig. 3. The graphs in Fig. 3 show both the 
upper and lower limit for the modulation loss as a function of RMS phase devi- 
ation 0. All information needed for analysis of the main carrier component is 
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RADIANS 


Figure 3-Upper and lower limits of modulation loss L m for a 
carrier phase modulated with RMS phase swing 0 

thus contained in these graphs. The use of the graphs in analysis work will be 
demonstrated in Chapter 5. 

It should once more be pointed out, that the results of this paragraph are 
valid only if the modulating frequencies are harmonically unrelated (see Equ. 
(3-3)). Modulation with harmonically related frequencies is treated in the next 
chapter. 


12 


4. Phase Modulation with Harmonically Related Frequencies 


If the modulating frequencies are harmonically related, we can find integers 
M n which are not all zero, so that 

= ± Mj ± M 2 co 2 ± ± M n a) n ± • • • ± M n = 0 (4-1) 

This implies that in calculating the carrier component from equ. (2-5) more 
terms then 

N 

A 0 = 7 ] ij ° 

1 

have to be included. For example, assume that the modulating function consists 
of three frequency components 

f (t) = <p 1 sin (cot + i/ij) + <p 2 sin(2c<Jt + \p 2 ) (4-2) 

+ cp 3 sin (3o)t +ifi 3 ) 

In this case the amplitude of the carrier component is 

\ - IJoCVJoWJoW -yo^Ji^JjCPi) 008 (20! -0 2 > 

+ 2J 1 (<p 1 )J 1 (<p 3 )J 2 ((p 2 )cos(0 1 - 20 2 + 0 3 ) 

- 2J 1 (cp 2 )J 1 (<p 3 )J 5 (<jp 1 ) cos (5\p x — 0 2 — 0 3 ) 

- j 2J 0 (9 3 ) Jj( 9 2 )J 2 ( 9 j) sin (20 2 - 0 2 ) 

+ J 210(93)12(92)14(01) sin ( 4 0 i " 20 2 ) 

- j 2 J 1 (9 1 )J 1 (9 2 )J 1 (<P 3 ) sin (0i + 02 - 0 3 > I (4 " 3) 


In the derivation, Bessel functions of higher order than 5 have been neglected. 
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Equ. (4-3) shows that with only three harmonically related modulating fre- 
quencies the expression for the carrier component already becomes a monstros- 
ity and a better method of analysis is desirable. If the modulating frequencies 
are harmonically related, they can be combined into a periodic function f (t) 


f (t) = f <t +T) 

The complex notation for the modulated carrier is 


y = e i(‘V+'/'o> e jf(t) 


(4-4) 


If f(t) is periodic, then ej f (*) is also periodic and the carrier component A 0 
can be found by complex Fourier analysis 


A 


o 



d t 


Because |e> f ( * > | < 1 the max value of A Q is 


(4-5) 


A < 1 (4-6) 

0 max — 

This max value is, for instance, obtained if f (t) is a narrow square pulse of 
width r as shown in Fig. 4, and we obtain 


1 im A Q = 1 (4-7) 

T — * 0 

independent of the RMS phase deviation 4 > . 

It has been proven in Appendix C that 

A Q > cos 0 (4-8) 

where 4 > is the RMS phase deviation, thus 


l 

T 


T 

J f(t) 1 2 dt=^ 2 
0 
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Figure 4— Narrow pulses for which Aq = Aq mQx = 1 if T-* 0 


The proof holds true if 


I f (t) | < 1.432 rr 


(4-9) 


and 


<t> < 1.432 77 


(4-10) 


Observing that cos 4 > becomes negative for 4 > > 77/2 we can combine Equ. (4-6) 
and (4-8) to 


and 


cos 4> < A n <1 for 4> < — 

_ 0 _ - 2 


A Q < 1 for <t> 1 2 


(4-12) 


It is also shown in Appendix C that A 0 min occurs if f (t) is an A-C s^nal with 
constant amplitude 


as illustrated in Fig. 5. 



(4-11) 
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Figure 5- Example of periodic time, function f(t) for which |f (t) | = 0 and 
*T 

f (t) dt = 0, so that A q = A q min = cos 0 


f 


The modulation loss L for the carrier component is thus 


and 


1 77 ""N 

10 log > L m > o for 4 > < — 

cos 2 4 > ^ 


> 


TT 


L > 0 for 0 > — 
m “ 2 > 


(4-13) 


Equ. (4-13) is shown graphically in Fig. 6. For comparison, also the modulation 
loss for the harmonically unrelated case is shown with dotted lines. We see that 
the upper and lower limits for the harmonically unrelated case are much closer 
together than for the harmonically related case. It is therefore of importance to 
distinguish between the harmonically related and unrelated case, especially for 
larger values of 0 . All A-C square waves, including pseudo random codes, and 
most amplitude limited signals satisfy Equ. (4-11) and must therefore be treated 
as harmonically related and the solid curves in Fig. 6 apply. An example where 
the modulating time function consists of both harmonically related and unrelated 
frequencies is treated in Chapter 5. 
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RADIANS 


Figure 6-Upper and lower limits of modulation loss L m fora carrier phase modulated 

with RMS phase swing <p 
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5. Applications 


Example 1 

As an example, let us analyze the case where a carrier is phase modulated 
with a time function consisting of a pseudo random code f 3 ( t) with amplitude g6, 
and two phase modulated subcarriers f 2 (t) and f 3 (t) with peak phase swings 4> 2 
and cf> 3 respectively. The total modulating time function f( t) is thus 

f(t) = f,(t) + f 2 (t) + f 3 (t) (5-1) 

We are interested in calculating the amplitude of the main carrier component coq 
for tracking purposes. Of interest are also the subcarrier components in the 
spectrum of the main carrier, i.e., a> 0 ± u> 2 and &> 0 ± o> 3 where a> 2 and co 3 are 
the angular frequencies of the subcarriers. We want to know the minimum of 
the main carrier components in order to insure receiver lock under worst con- 
ditions and the maximum of the subcarrier components in order to compute the 
risk of locking on to them instead of the main carrier. 

The phase modulated main carrier can be written 

V _ e i(“ot+tfo) e jf(t) 


or 


V = e H"0t +<i, 0) e jf 2 (t) e if 3 (t) 


(5-2) 


The Fourier expansion of the factor has a D-C component, whose 

magnitude A 01 is given by Equ. (4-8) 


A 01 = cos 4> x 


(5-3) 


because f j( t) is a square wave. The factors 
minimum D-C components A Q2 min and A 03 nin 


e and e J 

given by Equ. (3-7) 


have their 


A 02 min = 

A 03 „i„ = ^ 


(5-4) 


It is assumed that <£, < tt/2, $> 2 < 2.40 /V2, and ^> 3 < 2.40/^ so that Equ. (3-7) 
and (4-8) are applicable. If V given by Equ. (5-2) is expanded in a series, the 
coefficient for e i( " ot + ^ o) will be 


A 0 “ A 01 A 02 A 03 


(5-5) 
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and the minimum carrier component is thus 


A - A A A 

rt 0 min " n 01 min A 02 min rt 03 


( 5 - 6 ) 


or 


\ A Jo<^ ^2> 

The maximum modulation loss L m max corresponding to A 0 min is found from 
Fig. 6. With 


■= 0.6 rad 


^ = 0 3 =1.22 rad 


we obtain 



max 


max 


max 


= 1.7 dB 
= 3.6 dB 
= 3.6 dB 


and thus L =8.9 dB 

m max 


where L ml corresponds to A 01 etc. 
to Fig. 6, L, . = L. = 1.7 dB 

° 1 min 1 max 


L, . = 3.25 dB 

2 min 


L, . = 3.25 dB 

3 min 


L . = 8.2 dB 

m mm 


The minimum modulation loss is, according 
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The modulation loss for the main carrier will thus vary between 8.2 and 8.9 dB 
depending on the subcarrier modulation. Note that L is obtained for un- 
modulated subcarriers and L for modulated subcarriers. 

m min 

It has been assumed in the above derivation that the frequencies of fj(t) are 
harmonically related to each other and that the frequencies of f 2 (t) are not har- 
monically related to each other as well as the frequencies of f 3 (t). In addition, 
Equ. (5-5) is only true if none of the frequencies of fj(t) are harmonically re- 
lated to the frequency f 2 ( t ) or f 3 (t) and the frequencies of f 2 (t) are not har- 
monically related to the frequencies of f 3 (t) . 

If the frequencies are harmonically related the total RMS phase swing 4 > can 
be computed by 


* =}/'/’? (5-V 

and the carrier amplitude is then given by Equ. (4-8) 

A. . = cos 4 > (4-8) 

This is the absolute worst case. Using the same numbers as before we obtain 




= |/o.6* + I 1.22 


2 + il.22 ! 


= 1.363 rad 


and from Fig. 6 


L 

m max 


wo r s t 


case 


13.8 dB 


which is 4.9 dB more than the L m max computed under the assumption that the 
subcarrier and pseudo random code frequencies are harmonically unrelated. 

This example shows the importance of distinguishing between the harmonically 
related and unrelated case and of choosing suitable subcarrier frequencies. 

Let us also compute the subcarrier components in the spectrum of the modu- 
lated main carrier. 
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The amplitude A 2 of the frequency component co 0 ± co 2 is 


A 2 - A 01 A 12 A Q3 


(5-8) 


where 

✓ 


12 


J s Ccp,> 

W 


7 J 0^n) 


1 


is obtained from Equ. (A-ll) in Appendix A, and where cp n are the coefficients 
in the Fourier expansion of f 2 (t) . Again, Equ. (5-8) holds true only if the above 
assumption about the harmonic relation of the frequency components hold true. 

Because fj(t) is a square wave, A 01 is constant. A 12 max is obtained for 
<p x = /2 <p 2 and all other cp„ = 0, thus 

A ., 


03 max 


IS 


= e 

03 max 


A 


according to Equ. (3-7). Hence 


*3 


2 max 


= A 0J Jj(vT 4 > 2 ) e' 


The ratio of A, to A n is thus 

z max u mm 

1 1 

A 2 max Jl^ e 2 

A 0 m in " Jo(^ 


(5-9) 


(5-10) 
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and is obtained for subcarrier a> 2 unmodulated and subcarrier co 3 modulated. 
With the same numbers as before for <f> 1 , 0 2 , and 0 3 , we obtain 

A. 

-i- max = 0.794 = -2.0 db 
A. 

0 min 


The subcarrier components oj 0 ± a> 2 or a> 0 ± o> 3 may thus only be 2 db below 
the main carrier component o 0 . 


Example 2 

Assume that the RMS phase swing 0 is changed. How much does the mod- 
ulation loss L m change? 

For the harmonically related case we have (worst case) 

L m “ _2 ° lo S cos 4 > 

Differentiation yields 


dL 

— - = 20(log e) tan 0 db/rad (5-11) 

d 0 


For the harmonically unrelated cases we have (worst case) 


L m = -20 log J 0 (/2 0 ) 

and after differentiation 



20 (log e) 


S2 Jj(/2 0) 
J 0 (^2 0) 


(5-12) 


Again we consider the case where the modulating time function consists of 
a pseudo random code with RMS phase deviation 0 t and two subcarriers with 
RMS phase deviation 0 2 and 0 3 as in example 1. If the phase deviations have 
the same relative change 8 , we obtain 
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A0j = <f> x S 
A<^2 

A0 3 = 0 3 S 


(5-13) 


and from Equ. (5-11) and (5-12) to a first degree approximation 

AL ml = 8.69 ^ (tan <£j) S db 


^2 = 869 


/ 2 ^ J x (^2 0 2 ) 

J o (^2 0 2 ) 


S db 


(5-14) 


AL m3 = 8.69 


J ^/2 03 ) 


S db 


Assuming the worst case with all AL , , having the same sign, we thus obtain 

m 1*3 


for 


AL = AL . + AL , = AL , = 20S db 

m ml m2 m3 


cp 1 = 0.6 rad 


/2 0 2 = /2 4> 3 = 1.22 rad 


Example 3 


What is the spectrum of a phase modulated carrier, if f (t) is a symmetric 
square wave G(t) as shown in Fig. 7? 





T 


t 



1 

1 






Figure 7-Symmetric square wave 
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The modulated carrier is 


v = cos(^ 0 t + 0 O + G(t)) = Re{e i( " ot%) e’ 0 ^} (5-11) 

e j G ( t ) is a periodic function with period T and we can therefore expand 
into a complex Fourier series 


where 


and hence 


eiGCt) 



n=-oo 


c 

n 


i 

T 


j. 

T 



G(t) 







dt 


(5-12) 


(5-13) 




C 


sin nn cos <fi 
n 77 


+ 


1 


c-l) n 

nv 


sin 


Using complex notations we can thus write the modulated carrier 


V = g j ("o* +>/'o ) 


^COS + 



— sin zb e 

n77 



(5-14) 


where denotes summation over odd n only. 


This example demonstrates that it is easier in some cases to perform the 
complex Fourier analysis, Equ. (5-13), without expanding the modulating function 
into a Fourier series. Another example for which this is true, is a symmetric 
triangular modulating time function with peak amplitude $ and period T. The 
complex Fourier expansion yields 


V = e jC"ot +^o ) 


L 

n=— no 


2$ 

n 2 7 t 2 - 4<£ 2 




(5-15) 
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6. Generalizations 


So far the analysis has been concentrated on the component of carrier fre- 
quency. In this chapter an outline will be given how to analyze the other com- 
ponents in the spectrum. The harmonically unrelated and related cases will be 
treated separately as before. 


A. The Harmonically Unrelated Case 

Let us analyze the component of frequency co Q + co , where go. is given by 
Equ. (3-3). 


a; . = ± Mj a)j + M 2 a > 2 ± • • • ± M n co n ± • • • ± M N 

co n = angular frequency components of modulating function f(T) 

M n = integers, including 0 

The amplitude A. is obtained from Equ. (2-5) 

N 

~\ V-O t 6 ’ 1 ) 

nsl 

where cp n = the amplitude of component oj n 

If all M n = 0 then A. = A 0 . This case has been treated earlier. We therefore 
assume here that not all M n are zero. The lower limit, A . . , is zero, be- 

n i min 7 

cause if M ^ 0, then 



J„(<P) - 0 if <p - 0 


Thus 


A. . = 0 

i min 


(6-2) 


The maximum of A. may be calculated using Lagrange's multiplier with the 
side condition 



(6-3) 
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We obtain N equations 


3A. 

3 T + ^. = ° 


Performing the differentiation yields 


2K „ 

— t A -^ =0 
1 




(6-4) 


Together with (6-2), we have N + 1 equations, which is sufficient to eliminate 
the Lagrange's multiplier \ and solve for all cp n . Insertion in Equ. (6-1) will 
yield A i max . 


B. The Harmonically Related Case 

If the modulating time function f(t) is periodic with periodicity T, then a 
complex Fourier analysis of the modulated carrier is possible, resulting in 


where 


V = e^o^o) £c 




(6-4) 


(6-5) 


For some time functions such as square waves, triangular waves and sinusoidal 
waves, C n is easy to evaluate. This is demonstrated in Chapter 5, Example 3. 
Our problem is to find f(t) such that |C n | has a maximum or minimum. For 
n ^ Owe have 


because f(t) can be chosen so that all the power is in C 0 and therefore none in 

C . The upper limit, C has not been established so far. The problem is to 

n max 

find f ( t) , subject to the condition 

T 

f (t) 2 dt = 4> 2 



26 



such that |C n | given by Equ. (6-5) has a maximum and to calculate that 
maximum. 

The work is continuing along the above scheduled outlines (A) and (B) and 
more complete results will be published as soon as they are available. 

Another interesting generalization can be obtained from the fact that f (t) 
may not only be expanded into a Fourier series, but into any series consisting 
of orthogonal functions, provided, of course, that the set is complete. Let the 
orthogonal functions be g n (t) 


f (t) = 2] M 

0 


The modulated carrier is then 


V = e K«ot+*o> 


n,N 

[T e i<<Pn «n(t 

n.I 


)) 


(6-7) 


each factor e j(cPn Bn (t)) can be expanded in a complex Fourier series 


and 


e K<Pn8n(t>) 


max 00 



ms— oo 


e 


i 2nrr 

j T 


t 


(6-8) 


n.N 


V = 


J (" 0 t +'/'0 


’HE 


C e 

run 


j ^ nm7r t 

J X 


(6-9) 


If the g n (t)’s are sinusoidal functions, then the C nm 's are Bessel functions. If 
the g n (t)'s are square waves, the C nm 's are sine functions as shown in Chapter 5, 
Example 3. An analysis of phase modulation without involving Bessel functions 
is thus possible. 
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7. Conclusions 


It has been demonstrated that upper and lower limits of the spectral ampli- 
tudes of a phase modulated carrier can be calculated if the frequency spectrum 
of the modulating time function is not known and only the RMS phase deviation 
(= RMS of modulating time function) is known. Calculations can thus be carried 
out without introducing approximations such as substituting voice or other infor- 
mation by a single frequency. In addition, upper and lower limits are generally 
more important in circuit design and systems analysis than any value in between 
the limits, which exists only under very special conditions. The method is also 
suitable for stability and drift analysis, i.e., to analyze the change in the spec- 
trum of modulated carrier if the phase deviation changes. 
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APPENDIX A 


Phase Modulation with with Multiple Frequencies 

A carrier with angular frequency is phase modulated with N signals of 
angular frequency a ; and amplitude cp n> The instantaneous normalized ampli- 
tude v of the carrier is 


n=N 


v = cos|o> 0 t + </» 0 + Z n ^ ^ j 

or with complex notations 


n= 1 


/ n=N \ 

_ i +'A n )j 


V = e 


(A-l) 


(A-2) 


which also can be written 


where 


n 




y _ e i 0 t + ^0 ) 


n=N 


n 

n= 1 


e j<Pn sin(<u n t +</<„) 


denotes the product of the N factors e iCpnSin(< ° nt 


(A-3) 


Each of the factors is a periodic function in t and can therefore be expressed 
as a complex Fourier series 


e jCp sin(&)t +i£) 


z 


ms- oo 


c 

m 


e j m^t 


(A-4) 


Substituting o;t + 0 by 6 we obtain 

e jqp sin Q e jm(5-i />) 


(A-5) 


29 


where the complex constants C m are given by 


1 P + 7f 

_ 1 I j(Cpsin0-m# + 

= 2 tt 

J-TT 


mp) 


dd 


(A-6) 


From reference [7] we obtain 


1 f JC? si" 

2tt J 


6 - m0) 


de = JJ<9) 


(A-7) 


where J (cp) is the Bessel function of the first kind and order m. Hence 


C. = 


(A-8) 


and 


jm(<wt + 0) 


ei » .*»<«. ♦*> = jr J-(cp)e 

m=-co 

Substitution into Equ. (A-3) finally yields 

n=N m=i 

ttz: 

nsl m3— oo 

The carrier amplitude A 0 is obtained for m = 0, thus 


(A -9) 


n=N m=oo 

v,.'‘V^>TT 


(A-10) 


n»N 

*.-TTw 

n*l 


(A— 1 1) 


provided that the angular frequencies a> are not harmonically related. Under 
this condition the amplitude A nk for the frequency component a> 0 ± noj k is 


A 


nk 


J n K> 

J 0 C<P k ) 


n«N 

|7 

nsl 


(A-12) 


provided that naj fe cannot be obtained by adding or subtracting multiples of the 
other frequencies: 
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where M, etc are integers. 



APPENDIX B 


The problem is to find maxima and minima of 


n«N 

A, W <3 " 3) 

n=l 

where the independent variables <p n are subject to the condition 


nsN 

nssl 




Using the method of Lagranges multiplier, we obtain 




B A 

Bcp 


3A 0 B / 1 2 

B n + Bq>^2 Z,. “ 




n= 1 


n=N 


> 


0 2 =o 




where cp k and cp^ are any two of the variables cp n< Observing that 


r, ^ - - J ' M 


we obtain 


A 

V<P k ) ° 


+ A< P k 


= 0 




> 


-Ji(n) 

W 0 


+ = 0 


J 


(3-4) 


(B-l) 


(B-2) 
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and after elimination of A. and A„ 

Ji(9k) _ J 1 (y^) 

9 k Jo (9k) % Jo (9e) 

The obvious solution is 

9 k = 9e 

Because <p k and cp^ were arbitrarily chose, Equ. (B-4) holds for all <p n . A 0 has 
thus an extreme value for 

?» = To = ^ * (B-5) 


(B-3) 


tB-4t 


valid for all n . We still need to show that Equ. (B-3) has only one solution and 
we have to determine whether the extreme value of A 0 is a maximum, minimum 
or sadelpoint. By plotting 


J i (9) 
9 J 0 (9) 


(B-6) 


as is done in Fig. 8, we see that 


y = constant 

has only one solution <p <_ 5.1. By limiting <f> to 

0 < il* < JL2 n (B-6) 

"/2 " v /2 

Equ. (B-3) is thus limited to only one solution for all positive values of the 
integer N. 

In order to determine the type of extreme value of A Q , we study A 0 in the 
neighborhood of <p n = <p 0 by putting 

9 n = 9 0 + e n (B-7) 
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9 


where e n is an arbitrarily small quantity. From Equ. (3-4) we obtain 


n'N 

\ L (<P O + 2 <P 0 e « + e n)=^ 

n*l 

or 

n*N 

n*l n=l 

The expansion of J 0 (<p 0 + e n ) up to second order terms yields 

J<,(T„) = * <„> = J.(?o) - i - J 2 <*,>1 «J+- • 

2 \ T. 

Taking the logarithm of Equ. (3-3) yields 

n«N 

lnk 0 = J^ ^a,(f.)) 

n«l 

and substituting J 0 (? n ) from Equ. (B-9) 


¥ 


(B-8) 


(B-9) 


'in A 0 = N'fcn(J 0 (<p 0 ) ) + -!£*>) e n - if L<*'1 2+ . . . 

° \ Jo(9o ) ) J^oJoCP.) Ji(?o)/ 


For any I x | < | the logarithm 4n(l - x) can be expanded in a series 


(B-10) 


in (l - x) = - x - i x 2 . . . 

2 


and Equ. (B-10) can therefore be written, including up to second order terms 

n«N 


*n A 0 = N 4n(J 0 (<p 0 » V* e n 

JoW W 


Ji(?o) J 2 


Jo(To) 


n*=N 

X 

J n-1 
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* 


Substituting Ie n from Equ. (B-8) yields 


, = Nin(J 0 <?„)) i— [J,»o) 2 - Jo(fo)J 2 (?o)] r*2--- < B - ll > 

2J<,(T 2 ) 2 


In A 0 = 

Using the identities 


JoWo) = ^Ji»o) *5^ J lOPo> 

Jj(T 0 )=-Ji(%)-^- JiOPo) 

9 0 d 9o 


we obtain 


Ji(9 0 ) 2 - J 0 (9 0 ) J 2 (9 0 ) = 




(<Po> 


2 


The brackets [ ] in front of le* is thus positive for | <p 0 | > | . [ ] is also 
positive for | <p | <1, which can be shown in the following way: 


For cp 0 


= 0 we find 


[]- 


0. 


and -M ] = - Jo (<Po> J 2 (9 2 ) > 0 

d 9 0 9 0 

for | <p 0 | < | . As [] is an even function of <p 0 , we thus conclude that [ ] > 0 
for all cp Q . 

Hence, 


A q - e 


n= N 2 

“ C ?, £ n 

n = 1 


[J 0 (9 0 )] N 


(B-12) 


where C is a positive quantity. The extreme value is therefore a maximum. 
All cp n are equal at the maximum and we thus have 


A 


0 


max 


< 


lim 

N-00 



(B-13) 
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Expanding J 0 



into a series yields 


we also have 


with 



lim 

p— oo 



= e 


I (tl. *L + 

N \ 2 16N 



we obtain 



= lim e 
N— 00 


4 > 2 <£♦ 

2 + 16N 


e 


f. 

2 


(B-14) 


(B-15) 


The minimum value of A 0 has to be on the boundary because there was only one 
extreme value inside the boundary. The minimum value thus occurs for <p k = /2 0 
and all other <x> =0. Hence 

1 n 


for 


£ 

J o (/2 0) < A 0 < e" 2 
y/2 4> < 2.4048 


(B-16) 
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APPENDIX C 


O 


The problem is to find a periodic function f (t) such that A 0 has a mini- 
mum if A 0 is given by 





dt 


and where f (t) is subject to the conditions 


f 

Jo 


f (t) dt = 0 


and 


Using 


1 

T 



f (t) 2 dt = 4 > 2 


(C-l) 


(C-2) 


(C-3) 


e> x = cos x + j s in x 


Equ. (C-l) can be re-written 


1 

A o - y 


r T 

r T 

1 cos f (t) dt + j 

1 s in f (t) dt 


Jo 


and from 


we obtain 


|x + jy | - | x | 


J 

J 0 


*1 f 

T Jo 


cos f(t) dt >- cos f (t) dt 


(C-4) 
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is obtained. We consider all the f n to be independent variables and apply the 
method of Lagrange Multipliers 


and thus 



(C-7) 


-sin f k + 2\ f k = 0 


-sin fy + 2\ fy = 0 


(C-8) 


After elimination of A. 


sin f k sin fy 

fT 

k 


which has only one solution 


(C-9) 


ifj = i tt i 


(C-10) 


for 


|f k | < 1.43277 |fy| < 1.43277 

as can be seen from Fig. 10. 

From (C-6) follows that 


f n = ±4> (C-ll) 

Where half of the f n are positive and the other half are negative in order to 
satisfy Equ. (C-2). It thus follows that f(t) is an A-C signal with constant 
amplitude 0. For f(t) = ±0 we obtain 

1 r 

A Q = — I cos 0dt = cos <£ (C-12) 

Jq 
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or 


>1 f [l - JL<£2 + J_ 04 _ J_ «£6 
0 T J o L 2 ! 4 ! V 6 ! 


dt 


I f T \± .A 

T J o |_ 4 ! 6 ! 


0 2 ... 


e 2 dt 


(C-15) 


including up to e 2 terms. The e term vanishes due to Equ. (C-14). Thus 


A 0 ^cos</> + — I e 2 dt 


J r T 

r 

0 


where 


(C-16) 


K= 4? - 


4! 


A^ + ... ( — i > n (n+1) (n+2 -> + 

6! V ' 2(2n + 4)! 


(C-17) 


It is easily shown that 


1 /sinc/> 

K = ( - cos <p 

8 0 7 


2 V 0 


and 


K - 0 


for 0 ^ 1.4327T. Thus A 0 ^ cos 0 O for | 0 I - 1.432tt. 
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